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Abstract 

Massive stars end their life in an explosion event with kinetic energies of the order 10 51 erg. Immediately after the 
explosion has been launched, a region of low density and high entropy forms behind the ejecta which is continuously 
subject to neutrino heating. The neutrinos emitted from the remnant at the center, the protoneutron star (PNS), heat 
the material above the PNS surface. This heat is partly converted into kinetic energy and the material accelerates to an 
outflow that is known as the neutrino driven wind. For the first time, we simulate the collapse, bounce, explosion and 
the neutrino driven wind phases consistently over more than 20 seconds. Our numerical model is based on spherically 
symmetric general relativistic radiation hydrodynamics using spectral three flavor Boltzmann neutrino transport. In 
simulations where no explosions are obtained naturally, we model neutrino driven explosions for low and intermediate 
mass Fe-core progenitor stars by enhancing the charged current reaction rates. In the case of a special progenitor 
star, the 8.8 M Q O-Ne-Mg-core, the explosion in spherical symmetry was obtained without enhanced opacities. The 
post explosion evolution is in qualitative agreement with static steady-state and parametrized dynamic models of the 
neutrino driven wind. On the other hand, we find generally smaller neutrino luminosities and mean neutrino energies 
as well as a different evolutionary behavior of the neutrino luminosities and mean neutrino energies. The neutrino 
driven wind is proton-rich for more than 10 seconds and the contraction of the PNS differs from the assumptions 
made for the conditions at the inner boundary in previous neutrino driven wind studies. Despite the moderately large 
entropies of about 100 k B /baryon and the fast expansion timescale, the conditions found in our model are unlikely 
to favor r-process nucleosynthesis. The simulations are carried out until the neutrino driven wind settles down to 
a quasi-stationary state. After about 5 seconds post bounce, the peak temperature inside the PNS already starts to 
decrease due to the continued deleptonization. This moment determines the beginning of a cooling phase dominated 
by the emission of neutrinos. We discuss the physical conditions of the quasi-static PNS evolution and take the effects 
of deleptonization and mass accretion from early fallback into account. 
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1 Introduction 

Stars more massive than 8 M Q produce extended Fe-cores at the end of stellar evolution, which run into gravitational 
collapse due to pressure loss via photodisintegration and electron captures. The collapse halts at nuclear density, 
typically 2 - 4 x 10 14 g/cm 3 depending on the equation of state (EoS). However, the supersonically infalling material 
from the outer core continues to fall into the center. The core overshoots its equilibrium configuration and bounces 
back. A dynamic shock wave forms which propagates outwards and continuously loses energy due to the dissociation 
of heavy nuclei. As soon as the shock reaches the neutrinospheres, i.e. the neutrino energy and flavor dependent 
spheres of last scattering, additional electron captures emit a large amount of electron neutrinos. This burst of electron 
neutrinos, known as the deleptonization burst, carries away energy of several 10 53 erg/s on a timescale of 10 - 20 ms. 
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This energy loss turns the expanding shock into a standing accretion shock (SAS) already about 5 ms after bounce. 
Due to the large energy in the neutrino radiation field, neutrino heating between the neutrinospheres and the SAS has 



long been investigated as a possible source to revive the SAS and trigger neutrino driven explosions (Bethe and Wilson 
( [TM5| ), |Janka| ( [200T) , [Janka et al.| ( |2005] l, [Mezzacappa et al.| ( [2006l )). 

Up to now, neutrino driven explosions in spherical symmetry have only been obtained for the low mass 8.8 M Q O- 
Ne-Mg-core by Kitaura et al. ( 2006 ) and for low and intermediate mass progenitor stars by Sagert et al. ( 2009 ) assuming 
a hadron-quark phase transition during the early post bounce phase. On the other hand, multi-dimensional core collapse 
models with spectral neutrino transport have only recently become available and demonstrate the complexity of the 
underlying physical phenomena such as rotation and the development of fluid instabilities. These have been shown 
(see for example |Miller et al.~| ( [T993] l, |Herant et al~| (19941 ), |Burrows et IT] ( |T995] > and |Janka and Miiller |([T996ll) to 
increase the neutrino heating efficiency and help to understand aspherical explosions, see for example Bruenn et al 
(2006) and Marek and Janka (2009jl. For a review on axially-symmetric neutrino driven explosions, see also Janka 
|et al.||2008K 

The following dynamical evolution of the PNS and hence the properties of the neutrino spectra emitted is deter- 
mined by the mass accretion and the EoS. On a timescale of several seconds after the explosion has been launched, 
the region between the expanding explosion shock and the PNS at the center is subject to the formation of the neutrino 
driven wind as follows. Neutrinos continuously diffuse out of the hot PNS and heat the material on top of the PNS 
surface before they reach the neutrinospheres. We define the PNS surface to be the radius of the energy-integrated 
electron-neutrinosphere. The dominant neutrino heating contributions are given by the captures of electron-neutrinos 
and electron-antineutrinos on free nucleons. Matter is heated by neutrinos and the thermal energy is converted into ki- 
netic energy, which accelerates material on top of the PNS surface to positive velocities. This matter outflow is known 
as the neutrino driven wind. 

In this context, two particular studies are of special importance. The properties of the neutrino driven wind as 
described in Woosley et al. ( 1994) are based on the detailed radiation hydrodynamics simulation of a 20 M Fe-core 
progenitor applying the numerical model from Wilson and Mayle ( 1993). The simulation was carried until about 18 
seconds after bounce. Another state-of-the-art model of that time was the explosion of the O-Ne-Mg-core by Mayle 
and Wilson (1988]). Both investigations were milestones in the research of core collapse supernovae and are based 
on detailed neutrino input physics including neutrino transport, developed by J. R. Wilson. The results obtained, in 
particular the properties of the ejecta and the neutrino observables such as luminosities and energies, were considered 
the standard reference for more than 10 years. Neutrino driven wind studies used the results as parameters, where the 
conditions found indicated the possible site for the production of heavy elements via the r-process. In the simulations 



discussed in the present paper, we follow a similar approach as Wilson and Mayle ( 1993 i and Mayle and Wilson 
1988), where we apply the neutrino input physics based on Bruenn ( 1985[ ). Although we find general agreement with 



the previous work, in particular the explosion phase, we find significant differences in the properties of the neutrino 
driven wind. The entropies per baryon are smaller by a factor of 2-3 and the wind stays proton-rich for more than 10 
seconds for all our models. In addition, the neutrino luminosities and mean energies are generally smaller. The mean 
neutrino energies decrease with time, where they remain almost constant in the simulation of Woosl ey" et al.| (fl994). 
The largest difference arises in the decreasing difference between the mean electron neutrino and antineutrino energies 
found in our simulations, i.e. the neutrino spectra become more similar with respect to time. The difference in the 
neutrino spectra in Woosley et al. ( 1994 1 remains large and even increases with time. Using the results from Woosley 



et al.| ( [T994| ) as reference, Qian and Woosley ( 1996) analyzed the neutrino driven wind and formulated approximate 
analytical expressions for various properties of the neutrino driven wind, e.g. the neutrino heating rate, the electron 
fraction, the entropy per baryon and the mass outflow rate. 

Based on the static wind equations, the results obtained in parameter studies (see e.g. Duncan et al. ( 1986[ ), Hoffman 
et al. ( 1997[ ), [Thompson et al.|(|2001[ ) and Thompson and Burrows (2001 )) became known as static steady-state wind 
;re |\ 



models, where |Woosley and Baron| ( |1992| ), [Woosley et al.| ( |1994) , |Takahashi et al.| ( |1994| ) and |Witti et al.] ( |1994| ) 
described the neutrino driven wind in a radiation-hydrodynamics context. Of special importance for the neutrino driven 
wind investigations is the impact on the nucleosynthesis. Most interesting is the possibility to explain the production of 
heavy elements via the r-process due to the large entropies, the fast expansion timescales and the low electron fraction 
of Y e < 0.5 in the wind. Otsuki et al. ( 2000[ ) explored general relativistic effects on the neutrino driven wind and 
investigated the possible impact on the nucleosynthesis. Recently, Wanajo (2006a) and Wanajo (2006b) investigated 
the neutrino driven wind with respect to the r- and rp-processes. 



The possibility of supersonic wind velocities has been explored in most of the references. The supersonically 
expanding material in the wind collides with the much slower expanding and denser explosion ejecta. The material 
decelerates and a reverse shock forms which is known as the neutrino driven wind termination shock (first observed 
by Janka and Miiller| ( |1995[ ) and |Burrows et al" ( 1995| l). Recently, |Arcones et al. ( 2007| > examined the post bounce 
phase of core collapse supernovae of several massive progenitor stars. Their models were launched from massive 
progenitor stars that were previously evolved through the core collapse, bounce and early post bounce phases using 
sophisticated radiation hydrodynamics based on spectral neutrino transport in spherical symmetry. The simulations 
are then continued applying a simplified radiation hydrodynamics description (see Scheck et al. (2006l), assuming 
large luminosities to trigger neutrino driven explosions in spherical symmetry. The neutrino driven wind develops 
supersonic outflow and the wind termination shock appears in all their models. Like most of the present neutrino 
driven wind studies, an interior boundary was assumed instead of simulating the PNS interior for the PNS contraction 
and the diffusion of neutrinos out of the PNS. However, steady-state wind studies could not predict the important 
dynamical features from the presence of the neutrino driven wind termination shock, especially the deceleration of the 
wind material and the consequent entropy as well as density and temperature increase during the deceleration. In this 
respect, the investigation from |Arcones et aL] ( |2007| ) was a milestone in modeling the neutrino driven wind consistently. 
On the other hand, they were focusing on parameters (luminosities and mean neutrino energies) in agreement with the 
simulations of Bethe and Wilson (19851) and Woosley et al. ( 1994 1. Although the dynamics is in general agreement, 
several properties of the neutrino driven wind as well as the neutrino spectra emitted differ significantly. 

The present paper follows a different approach. We simulate consistently the dynamical evolution of the collapse, 
bounce and post bounce phases until the neutrino driven wind phase for more than 20 seconds. The simulations are 
launched from the 8.8 Mq O-Ne-Mg-core from Nomoto (1983,1984,1987) and the 10.8 and 18 M Fe-core progenitors 
fromfWoosley et al. (2002). Our numerical model is based on general relativistic radiation hydrodynamics with spec- 
tral three-flavor Boltzmann neutrino transport in spherical symmetry. The explosion mechanism of massive Fe-core 
progenitors is an active subject of research. To model neutrino driven explosions for such progenitors in spherical 
symmetry, we enhance the electronic charged current reaction rates artificially which increases the neutrino energy 
deposition and revives the SAS. The mechanism including the tuned neutrino reaction rates will be further discussed in 
§2 and §3. Such explosion models were investigated with respect to the nucleosynthesis by Frohlich et al. (2006 a-c). 
Here, we report on the formation of the neutrino driven wind and the possibility of the wind developing supersonic 
velocities and hence the formation of the wind termination shock. In addition, we will also illustrate the explosion and 
the neutrino driven wind phases of the O-Ne-Mg core, where an explosion is found in spherical symmetry even without 
artificially enhanced neutrino opacities. The results are in qualitative agreement with those of Mayle and Wilson ( 1988 ) 
and |Kitaura et al. (2006), who used a different EoS. 

We find that for low progenitor masses, the neutrino driven wind termination shock will develop, using the tuned 
neutrino reaction rates. When the neutrino reaction rates are switched back to the standard opacities given in|Bruenn 



(1985]), the neutrino driven wind develops only a subsonic matter outflow. For intermediate progenitor masses, the 
neutrino driven wind remains subsonic even with the artificially enhanced neutrino emission and absorption rates. 
Since the neutrino driven wind depends sensitively on the emitted neutrino spectra at the neutrinospheres, we believe 
accurate neutrino transport and general relativity in the presence of strong gravitational fields are essential in order to 
describe the dynamical evolution. Furthermore, the accurate modeling of the electron fraction in the wind is essential 
for nucleosynthesis calculations, which can only be obtained solving the neutrino transport equation. In addition, it is 
beyond the present computational capabilities to carry multi-dimensional simulations with neutrino transport to several 
seconds after bounce. Hence, our investigations are performed in spherical symmetry where we simulate the entire 
PNS interior rather than approximating an interior boundary. We find significant discrepancies in comparison with the 
assumptions made in previous wind studies. Material is found to be proton rich for more than 10 seconds, where most 
wind models assume luminosities and mean neutrino energies such that the neutrino driven wind becomes neutron-rich. 
We question the validity of the approximations made in such wind studies. We believe that an accurate and consistent 
modeling of the physical conditions in the wind is essential, especially in order to be able to draw conclusions with 
respect to nucleosynthesis. 

The paper is organized as follows. In §2, we will present our spherically symmetric core collapse model. §3 is 
devoted to the explosion phase of neutrino driven explosions in spherical symmetry. We examine the 8.8 M Q progenitor 
model from Nomoto (1983,1984,1987) using the standard neutrino opacities and the 10.8 and 18 Mq progenitor models 
from Woosley et al. (2002]) using artificially enhanced neutrino reaction rates. In §4 we discuss the conditions for the 



formation of the neutrino driven wind and the possibility for the wind to develop supersonic velocities. We discuss in 
§5 the electron fraction approximation used in the literature. Since a generally neutron-rich neutrino driven wind is 
found in many previous and present wind studies, we illustrate the differences and investigate why we find a generally 
proton-rich wind. §6 is dedicated to the long term post bounce evolution for more than 20 seconds. In §7 we discuss 
the results and emphasize the main differences of the present investigation to previous wind studies. Finally we close 
with a summary in §8. 



2 The model 



Our core collapse model Agile-Boltztran is based on general relativistic radiation hydrodynamics in spherical symme- 
try, using three-flavor (anti)neutrino Boltzmann transport. For details see Mezzacappa and Bruenn (1993 a-c), Mezza- 
cappa and Messer ( 1999), Liebendorfer et al. (2001a,b) and Liebendorfe r et al.| ( |2004| l. For this study we include the 
neutrino input physics based on Bruenn ( 1985| l. The charged current reactions considered 



e + p ^ n + v e , 

e + + n ^ p + v e , 
e~ +{A,Z) ^ <A,Z-1> 



(1) 
(2) 
(3) 



are electron and positron captures at free nucleons as well as electron captures at nuclei. The nuclei are characterized 
by an average atomic mass and charge (A, Z). In addition, the standard scattering reactions considered are iso-energetic 
neutrino nucleon (N e {n,p}) and nuclei (N = (A,Z)) scattering, 



V + N ■ 



■N, 



where v e {v e , v^/ T } (equivalent for antineutrinos v), and neutrino electron/positron scattering 

v + e ± ^ v + e ± . 

The classical neutrino pair process is electron-positron annihilation, 

e~ + e~ ^ v + v. 

The standard neutrino energy E dependent emissivity j(E) and absorptivity x(E) for the charged current reactions as 
well as the scattering and pair-reaction kernels are given in Bruenn ( 1985 1 based on Yueh and Buchler ( 1976} and 
Schinder and Shapiro ( 1982[ l. The additional pair-process nucleon-nucleon-Bremsstrahlung, 



N + N^N + N + v + v, 



has been implemented into our model according to Thompson and Burrows (2001 1 and is also taken into account. The 
emission of ///r-neutrino pairs via the annihilation of trapped electron-neutrino pairs, 



v« + v. 



+ V 



V/T, 



as well as contributions from nucleon-recoil and weak magnetism as studied in Horowitz ( 2002[ l are investigated in 
Fischer et al. ( 2009[ l and are not taken into account in the present study of the neutrino driven wind. 



2.1 Recent improvements of the adaptive grid 

Long-term simulations of the supernova post bounce phase with Agile-Boltztran lead to a very large contrast of densi- 
ties, reaching from ~ 10 15 g/cm 3 at the center of the protoneutron star (PNS) to densities of order g/cm 3 and smaller in 



the outer layers. The version of Agile described in Liebendorfer et al. (2002i is not able to resolve such large density 
contrasts. If the enclosed mass a is large and the density in one zone very small, then the evaluation of the mass 
contained in the zone according to Eq. (39) in Liebendorfe r et al.] ( |2002| l, 



da i+ i = a M - a h 



is subject to large cancellation so that truncation errors inhibit the convergence of the Newton-Raphson scheme in the 
implicitly finite differenced time step. However, the problem can be avoided by a simple modification that was first 
explored in Fisker (2004). The state vector of Agile-Boltztran is given by the following set of quantities 



y — (a , r ,u ,m ,p ,T , Y e ) , 



(4) 



with enclosed baron mass a, radius r, velocity u, gravitational mass m, baryon density p, temperature T, electron 
fraction Y e . In the improved version, the state vector at time f is based on zone masses, da" , , where the enclosed 
mass 



/-l 

a'; - ) da" 



! 



1 



becomes the derived quantity. 

The form of the generic equation (30) in Liebendorfer et al. ( 2002 1 applies to the continuity equation, the momen- 
tum equation and the energy equa tion. If we define < 5, - a n+1 — a" as the difference of the enclosed mass a, between 
time f"and f" +1 , Eq. (30) in Liebendorfer et al. (2002i becomes 



'I (da" + 6 M -6i)-tf 

2 \ 2 I H 



da" . 

i+4 i+i 



dt 



where the relative velocity between the fluid and the grid in the advection term y is defined by 



dt 



A 

dt' 



(5) 



(6) 



The cancellation of large numbers during the Newton-Raphson iterations of the implicit time step is avoided if the time 
shifts <S, are chosen as the unknowns in the state vector when Eqs. |5]l and |6]i are solved. The vector of zone masses is 
then updated between the implicit time steps by 

da" + j - da" , + 6m - 6i. 

This leads to satisfactory convergence of the Newton-Raphson iterations even in the presence of large density contrasts. 



2.2 The equation of state 



For the present investigation of the neutrino driven wind, the baryon EoS from Shen et al. ( 1998 ) for hot and dense nu- 
clear matter has been implemented for matter in nuclear statistical equilibrium (NSE). For temperatures below T = 0.5 
MeV where NSE does not apply, the baryon EoS combines an ideal gas approximation for a distribution of nuclei, 
based on Tim mes and Arnett| ( |1999| l (including ion-ion-correlations), and a nuclear reaction network using the com- 
position given by the progenitor model. Details of the reaction network are given in Thielemann et al. ( 2004 1 and 
references therein. The nuclear abundances are included in the state vector of Agile-Boltztran as follows 



y — (a , r , u , m , p , T , Y e , Y\ , . . . , Y^) ■ 



(7) 



For all these quantities, including the nuclear abundances Y\, Y^, the corresponding advection equations are solved 



as described in |Liebendorfer et al.| l |2002) §3, but with an improved second order accurate total variation diminishing 
advection scheme. The nuclear reaction network is used in an operator-split manner in order to calculate the abundance 
changes due to the source terms which in turn depend on employed reaction rates. 

Due to computational limitations, we restrict ourselves to N = 19. We consider the free nucleons and the 14 sym- 
metric nuclei, from 4 He to 56 Ni. In order to model matter with Y e + 0.5 to some extent, we additionally include 53 Fe, 
54 Fe and 56 Fe. The network calculates the composition dynamically from the progenitor stage until the simulations are 
stopped. It is used for an accurate internal energy evolution. In addition, we can approximately reflect the composition 
of the PNS surface for more than 20 seconds after bounce, where nuclei that have been previously in NSE are freezing 
out of NSE as the temperature drops rapidly below 0.5 MeV already about 1 second after bounce and reach below 0.01 



MeV at about 10 seconds post bounce. In previous studies the simplification of an ideal gas of Si-nuclei was used for 
matter which is not in NSE. This leads to an increasing inaccurate internal energy evolution after 500 ms post bounce 
when the explosion shock reaches the Si-layer of the progenitor and simplifications could not be extended beyond 1 
second post bounce time. The implementation of the nuclear reaction network now makes it possible to include more 
mass (up to and including a large fraction of the He-layer, depending on the progenitor model) into the physical domain 
and follow the dynamical evolution of the explosion by one order of magnitude longer. 

The baryon EoSs for NSE and for non-NSE are coupled to an electron-positron EoS (including photons), developed 



by |Timmes and Arnett| (|1999). 



2.3 Enhanced neutrino emissivity and opacity 

By our choice of a spherically symmetric approach, we implement the explosion mechanism of massive Fe-core pro- 
genitor stars artificially to trigger neutrino driven explosions during the post bounce evolution after the deleptonization 
burst has been launched. Neutrino heating between the neutrinospheres and the SAS transfers energy from the radia- 
tion field into the fluid. A part of this energy is converted into thermal energy which revives the SAS and launches the 
explosion. The revival of the SAS and hence the neutrino driven explosions take place on a timescale of several 100 
ms. 

During the post bounce evolution, heavy nuclei continue to fall onto the SAS and dissociate into free nucleons. 
These free nucleons accrete onto the PNS surface. Hence the dominant neutrino heating contributions behind the SAS 
are due to the electronic charged current reactions Eqs. ([T]) and (|2]). To trigger explosions in spherically symmetric 
core collapse simulations of massive Fe-core progenitors, we enhance the emissivity j and absorptivity x by a certain 
factor (typically 5 - 7) in the region between the SAS and the neutrinospheres. This corresponds to matter with 
entropies above 6 kfi/baryon and baryon densities below 10 10 g/cm 3 . The entropies ahead of the shock are lower 
and the central densities of the PNS are higher, such that the artificial heating only applies to the region between the 
neutrinospheres and the SAS. The artificially enhanced reaction rates do not change the neutrino luminosities and 
mean neutrino energies significantly for the electron-neutrinos and electron-antineutrinos. Furthermore, /^-equilibrium 
is fulfilled since the reverse reaction rates are obtained via the detailed balance. However, the timescale for weak- 
equilibrium to be established is reduced and hence the electron fraction changes on a shorter timescale. In combination 
with the increased neutrino energy deposition, this leads to a deviation of the thermodynamic variables in comparison 



to simulations using the standard opacities given in Bruenn ( 1985 1, which will be further discussed in §7. The weak 
neutrino driven explosions obtained have explosion energies of 6.5 x 10 50 erg and 2 x 10 50 erg for the 10.8 and the 18 
M Q progenitor model respectively. 

2.4 Explosion energy and mass cut 

The explosion energy estimate is a quantity which changes during the dynamical evolution of the system. It is given 
by the total specific energy of the fluid in the laboratory frame 

£s P ecific(f, a) = Ye + ^-j (y - ™) . (8) 

which in turn is the sum of the specific internal energy e [j] the specific kinetic energy given by the fluid velocity 
u — dr/adt squared and th e specific gravitational binding energy m/r with gravitational mass m and radius r (see 



Liebendorfer et al 



( 2001b I). Y(t,a) = Vl - 2m/r + u 1 and a(t, a) are the metric functions in a non-stationary and 
spherically symmetric spacetime with coordinate time f, baryon mass a and the two angular coordinates (0, <p) describ- 
ing a 2-sphere of radius r(t, a) [see Mis ner and Sharp| ( |1964[ )]. The explosion is determined by the energy of the ejecta. 
Integrating E spec j& c (t, a) with respect to the enclosed baryon mass starting from the progenitor surface M toward the 
center 

£totai(fo> «o) = _ I E specl f lc (to,a)da, (9) 



The baryon contribution to the internal energy is composed of a thermal and nuclear part, i.e. e = ethermal + ^nuclear- In NSE, e is given implicitly 
via the EoS of hot and dense nuclear matter. In non-NSE, e nuc iear is the binding energy of the nuclear composition used in the reaction network. 



gives the total mass-integrated energy, at a given time to outside a given mass «o- The expression (|9]l is negative during 
the collapse, bounce and the early post bounce phases because the progenitor and central Fe-core are gravitationally 
bound. At some time after bounce, expression (|9]l becomes positive in the region between the shock and the neutri- 
nospheres. It stays negative at large distances and close to the deep gravitational potential of the PNS, because the 
outer layers of the progenitor and the PNS continue to be gravitationally bound. While the emission of neutrinos cools 
the PNS interior, neutrino absorption deposits energy on a timescale Seating of the order 100 ms into the fluid near the 
neutrinospheres. This increases the specific internal energy which matches at later (~ 500 ms) post bounce times the 
gravitational binding energy at a certain distance toward the center. 
On a suggestion by S. Bruenn, we consider the mass cut as follows 

"cut = a (max \E total (t » T heating , a))) . (10) 

The material outside of a cut is gravitationally unbound and will be ejected while the enclosed material will accrete onto 
the central PNS. The explosion energy estimate is defined as the total mass-integrated energy of the layers outside the 
mass cut 

^expl = £total(f ^ Theating > Acut) > (11) 

at late times (t » Theating) a ft er the explosion has been launched. It becomes clear from the above expressions that 
the explosion energy estimate is sensitively determined by the balance of internal and kinetic energies to gravitational 
binding energy. 

From the time post bounce when the shock reaches low enough densities and temperatures such that neutrinos 
decouple from matter completely, neutrino heating and cooling does not affect the explosion energy estimate anymore. 
The additional energy deposition from the neutrino driven wind, which will be discussed further below, might affect 
the explosion estimate at later times. We will illustrate in particular the effect of the formation of a supersonic neutrino 
driven wind and the wind termination shock to the explosion energy estimate. Only after the neutrino driven wind 
disappears, the final value of the explosion energy can be obtained. 



2.5 The neutrino observables 

The neutrino radiation hydrodynamics equations are a coupled system which combines the evolution of the matter 



properties given by the state vector y and the radiation field, as documented in Liebendorfer et al. ( 2004 1 and references 
therein. The neutrino radiation field is taken into account via the phase-space distribution function f v (t, a,fj.,E) for each 
neutrino flavor v = (v e , v e , v^ T , v^/j). In spherical symmetry, it depends on the time f, the enclosed baryon mass a as 
well as on the neutrino energy E and the cosine of the propagation angle /i = cos(#). The evolution of the distribution 
function is modeled via solving the Boltztran transport equation. It determines the phase-space derivative of the specific 
distribution function F v - f v /p, i.e. the distribution function divided by the matter density p, in a co-moving frame 
(see for example Eq.(8) in [Liebendorf er et al.| ([T005 )) and due to neutrino-matter interactions such as emission and 
absorption as well as scattering and pair reactions. 

In order to compare simulation results, neutrino observables can be defined. Commonly used is the neutrino 
number-luminosities, which is given by the phase-space integration of the neutrino distribution function as follows 

2 2nc f +1 r°° 2 
L„ (Uo, h],a) — Anr p— —? I da I E dE F v (t,a,u,E), 
(hey J_i Jo 

which is the number of neutrinos^] of energy E passing through the mass coordinate a for a given time-interval [?o,fi] 
taken in a co-moving frame at position r(t, a). Equivalently, the energy-luminosity can be defined as follows 

L e ([t ,ti],a) = Ly([t ,h],d) 

= Anr 2 p r- udp. E dE F v (t,a,p,E), 

(hey J-i Jo 



The integration with respect to /i is performed separately for in- and out-ward direction, according to the transport coefficients. 



for each neutrino flavor, i.e. (v e , v e , v„/ T , v^). Additionally useful quantities are the mean neutrino and root-mean- 
squared (rms) neutrino energies, defined as follows 



(E v (t,a)) 



(E v (t, a)) n 



J^ 1 dfi E 3 dE F v (t, a, 11, E) 
jTj 1 dfi £° E 2 dE F v (t, a, 11, E) ' 



JT, dfi E 4 dE F v {t, a,n,E) 
\ dfx J^°° E 2 dE F v (t, a,n, E) 



We will use these observables, i.e. the energy-luminosities and the mean and root-mean-squared energies, to illustrate 
the dynamical evolution of the radiation field as well as for comparisons with previous studies. 



2.6 The electron fraction 

The proton-to-baryon ratio is essential for the composition of the ejecta, which is obtained via detailed post processing 
nucleosynthesis calculations. In the absence of muons or tauons, the proton-to-baryon ratio is given by the electron 
fraction as follows 

Y e = Y e --Y e+ , (12) 

which defines the total number of charges per baryon. The change of the electron fraction is given by the balance of 
emitted and absorbed electrons (positrons) and electron-neutrinos (antineutrinos) at free nucleons and nuclei. Weak- 
equilibrium is achieved if 

He- + Hp = fin + fX Ye , (13) 



(14) 



where ju, are the chemical potentials for electron and positron (ju e ±), proton (fi p ), neutron (ji„) and electron-neutrino 
(p-nu e ) and electron-antineutrino QiyJ. The time-derivative of the electron fraction, Y e , is given by the phase-space 
integration of the emissivities j v and the absorptivities Xv f° r electron-neutrinos and electron-antineutrinos as follows 



Ye = - 



2n nig c 
(he) 3 ~ 



E dE x 



X (O'v, -Xvefve) ~ O'v- ~ Xv e fvJ) 



(15) 



where m B is the baryon mass, p is the matter density and^ = j + X- The emissivities j v (E) and absorptivities Xv(E) 
depend on the neutrino energy. They are the reaction rates for the electronic charged current reactions, which are 
calculated following Bruenn (1985) and depend on temperature and density. Eq. (15i is found by combining the 
equation of lepton number conservation with the phase-space integrated Boltzmann transport equation (see Mezzacappa 
& Bruenn (1993 a)). To follow the dynamical evolution of the electron fraction via Eq. ( |15) , knowledge of the neutrino 
distribution functions is required for which neutrino transport is necessary. 



3 Explosions in spherical symmetry 

Progenitor stars more massive than 9 M Q develop extended Fe-cores at the end of stellar evolution. The explosion 
mechanism of such Fe-core progenitors is an active subject of research. In the following section, we will investigate 
the neutrino driven explosions of the 10.8 and 18 M Q Fe-core progenitors from|Woosley et aT (2002 1 in spherical 



symmetry by enhancing the electronic charged current reaction rates artificially. Further below, we will investigate the 
explosion phase of the 8.8 M Q O-Ne-Mg-core from Nomoto (1983,1984,1987), where the explosion is obtained using 
the standard neutrino opacities as introduced in § 2. 



3.1 Neutrino driven explosions of Fe-core progenitors 

The dynamical behavior of the explosion energy estimate and the shock position are illustrated in Figs.[T](a) and (b) 
respectively with respect to time after bounce. The figures illustrate the explosion phase and long time evolution up to 
10 seconds after bounce. After achieving a convergent value between 600 ms and 2 seconds post bounce of 4.5 x 10 50 
erg, the explosion energy estimate is lifted slightly to about 6.5 x 10 50 erg. This effect coincides with the additional 
mass outflow obtained in the neutrino driven wind and the appearance of the reverse shock, which will be discussed 
further below. In simulations with a less intense (subsonic) neutrino driven wind, this effect is much weaker and the 
explosion energy can be obtained already after about 1 second post bounce. 




Figure 1: Explosion energy estimate and shock position with respect to time after bounce for the 10.8 M progenitor 
model from Woosley et al. ( 2002 1. In addition, graph (b) illustrates the position of the He-layer (dashed line) and the 
O-layer (dash-dotted fine). 



The neutrino luminosities and the mean as well as rms neutrino energies are shown in Fig|2]for the 10.8 (middle 
panel) and the 18 (right panel) M Q progenitor model with respect to time after bounce. Note that the more compact 
PNS of the 18 M Q progenitor model results in generally larger neutrino luminosities. The oscillating shock position 
and the consequent contracting and expanding neutrinospheres during the neutrino heating phase of the 10.8 and 18 
M Q progenitor models on a timescale of several 100 ms are reflected in the electron-flavor neutrino luminosities, 
which correspondingly increase and decrease respectively. During the heating phase, the mean neutrino energies of 
the electron-(anti)neutrinos increase from about 8 (10) MeV after bounce to about 12 (14) MeV until the explosion is 
launched. The mean neutrino energy of the ^/r-neutrinos remains constant at about 18 MeV during the heating phase. 
The mean neutrino energies are generally smaller than the rms-energies but follow the same behavior. The explosions 
for both models are launched after about 350 ms post bounce, which defines the neutrino heating timescale for the 
energy deposition in the gain region to revive the SAS. Matter is accelerated to positive velocities and the SAS turns 
into the dynamic explosion shock. The resulting neutrino spectra from these artificially induced explosions in spherical 
symmetry are in general agreement with the neutrino spectra from axially-symmetric neutrino driven core collapse 
supernova models that explode without artificially modified reaction rates (see Marek and Janka (2009)). The explosion 
shock continuously propagates through the remaining domain of the progenitor star. After the explosions have been 
launched, the electron flavor neutrino luminosities decay exponentially. Furthermore the jumps in the neutrino energies 



after 350 ms post bounce for the 10.8 and 18 M progenitor models are due to the shock propagation over the position 
of 500 km, where the observables are measured in a co-moving reference frame. 
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Figure 2: Neutrino luminosities and energies with respect to time after bounce for the 8.8 M Q O-Ne-Mg-core progen- 
itor model from Nomoto (1983,1984,1987) (left panel) and the 10.8 M Q and 18 M Q Fe-core progenitor models from 



Woosley et al. (20021 (middle and right panels respectively), measured in a co-moving frame at 500 km distance. 



3.2 The O-Ne-Mg-core 

A special star is the 8.8 M Q progenitor model from Nomoto (1983,1984,1987). The central thermodynamic conditions 
at the end of stellar evolution are such that only a tiny fraction of about 0. 15 M Q of Fe-group nuclei are produced, where 
nuclear statistical equilibrium (NSE) applies (see Fig. [3] (a) top panel). Instead, the central composition is dominated 
by 16 0, 20 Ne and 24 Mg nuclei. Because temperature and density increase during the collapse, these nuclei are burned 
into Fe-group nuclei and the NSE regime increases (see Fig.|5]middle panel). The core continues to deleptonize, which 
can be identified at the decreasing Y e in Fig. [3] We use our nuclear reaction network as described in §2.2 to calculate 
the dynamically changing composition, based on the abundances provided by the progenitor model. The size of the 



bouncing core of M cole 0.65 M is significantly larger in comparison with the previous studies by Kitaura et al. 
(2006) and Liebendorfer (2004i, illustrated in Fig. 4(a) at different velocity profiles before and at bounce. This is 
because we do not take the improved electron capture rates from Hix et al. ( 2003 \ and Langanke et al. ( 2003 1 into 
account, which are based on the capture of electrons on the distribution of heavy nuclei. It results in a lower central 
electron fraction at bounce and a consequently more compact bouncing core of about ^0.1 M , in comparison to the 
standard rate s given in | Bruenn ( 1985[ l. The remaining difference is most likely due to the large asymmetry energy of 
the EoS from Shen et al. ( 1998 1 applied to the present study. 

This progenitor is not only a special case due to the incomplete nuclear burning at the end of stellar evolution but 
also due to the steep density gradient which separates the dense core from the He- and H-envelope at 1.376 M , see 
Fig.[3](c). There, the density drops over 13 orders of magnitude which makes it difficult to handle numerically. 



10 



.2 10~ 



10"' 



10 





10 



o 10" 



a io"- 



10~ 







1 

uniform , , 





nuclear !!! 




matter "} , 
i i't - 


1 1 


it 

:l 'i 

1 1 

1 1 
1 i 


i 

j | 
i i 
i 1 
1 1 

j 



0.5 



0.4 



E 0.3 



0.2 



0.2 0.4 0.6 0.8 1 1.2 1.4 
Baryon Mass, a [solar masses] 



0.1 



0.5 



c 0.4 



S5 0.3 



0.2 



ru 0.1 



0.2 0.4 0.6 0.8 1 1.2 1.4 
Baryon Mass, a [solar masses] 



0.2 0.4 0.6 0.8 1 1.2 1.4 
Baryon Mass, a [solar masses] 



0.2 0.4 0.6 0.8 1 1.2 1.4 
Baryon Mass, a [solar masses] 



0.5 



c 0.4 



2 0.3 



0.2 



Qj 0.1 



14 
12 

10 : 



fe- 
rn 

CO 



f; 14 



0.2 0.4 0.6 0.8 1 1.2 1.4 
Baryon Mass, a [solar masses] 



0.2 0.4 0.6 0.8 1 1.2 1.4 
Baryon Mass, a [solar masses] 



-progenitor 



0.2 0.4 0.6 0.8 1 1.2 1.4 
Baryon Mass, a [solar masses] 

24 ms before bouncel 




0.2 0.4 0.6 0.8 1 1.2 1.4 
Baryon Mass, a [solar masses] 




0.2 0.4 0.6 0.8 1 1.2 1.4 
Baryon Mass, a [solar masses] 



Figure 3: Evolution of the 8.8 M progenitor model from Nomoto (1983,1984,1987) during the core collapse phase 
(top: progenitor configuration, middle: 24 ms before bounce, bottom: at bounce). The composition (left panels) are 
as follows: C+O (thin solid line), Ne+Mg (dashed line), Fe+Ni (dash-dotted line), He (dotted line), Si+S (thick solid 
line). 



The low density of the mass outside the O-Ne-Mg-core makes it possible to obtain the explosion in spherical 
symmetry supported via neutrino heating. Neutrino cooling in the region of dissociated nuclear matter causes the 
expanding shock front to turn into the SAS with no significant matter outflow, v^-cooling dominates over v e -heating by 
one order of magnitude. Only at the dissociation line of infalling heavy nuclei, the neutrino energy deposition drives 
the SAS slowly to larger radii, for illustration see the heating(cooling) rates and velocity profile in Fig.[5](left panel) at 
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Figure 4: Radial velocity profiles with respect to the baryon mass (a) and with respect to the radius (b) for the 8.8 M G 
progenitor model from Nomoto (1983,1984,1987) 



20 ms post bounce. However, the cooling of v e still contributes to a large amount at 25 ms post bounce over the heating 
of v e and v A ,/ T in Fig. [5] (middle panel) behind the SAS. Only directly at the shock a small net-heating rate remains. 
Hence the influence of the neutrinos to the explosion itself is of minor importance. More important is the region of C- 
O-burning which produces Ne and Mg. The hydrodynamic feedback to this thermodynamic transition can be identified 
already during the collapse phase of the progenitor core at the velocity profiles in Fig. 4(a) at about 1.35 - 1.374 M . 
As material is shock heated post bounce, the transition layer where Ne and Mg nuclei are burned into NSE propagates 
together with the expanding shock wave outwards. In other words, the Ne-Mg-layer of the progenitor is converted 
directly into NSE. Furthermore, the transition (discontinuity) from C-O-burning is falling quickly towards the SAS. It 
was found to be at about 350 km at 20 ms post bounce and at about 200 km at 25 ms post bounce, illustrated at the 
velocity profiles (bottom) in Fig. [5] (left-right panels). At about 30 ms post bounce, the entire Ne-Mg-layer is converted 
into NSE due to the temperature increase obtained via shock heating. Hence, C and O nuclei are burned directly into 
NSE. 

In contrast to more massive Fe-core progenitors where O-burning produced an extended Si-S-layer, the amount of 
28 Si and 32 S is less then 1% at the end of nuclear burning for the O-Ne-Mg-core discussed here (see Fig. [3](a) middle 
panel). This small fraction of Si and S is already converted into NSE during the initial collapse phase, due to the rapid 
density and temperature increase of the contracting core. Hence, C- and O-nuclei are burned directly into NSE during 
the post bounce evolution. This sharp transition is related to a jump in the density and the thermodynamic variables. As 
the SAS propagates over this transition along the decreasing density, the shock accelerates to positive velocities (see 
Fig. |5]right panel). The consequent explosion is hence driven due to the shock propagation over the infalling transition 
between two different thermonuclear regimes rather then by pure neutrino heating, illustrated at the velocity profiles 
in Fig. 4(b) Although Kitaur a~et al.| ( |2006| approximated nuclear reactions during the evolution of the O-Ne-Mg-core 
progenitor, the results of their explosion dynamics are in qualitative agreement with our findings. 

The more massive Fe-core progenitors show the same thermo- and hydrodynamic features as discussed above for 
the O-Ne-Mg-core due to the transition between different thermonuclear regimes. However, the differences are smaller 
because C-O-burning produces an extended layer composed of 28 Si and 32 S. The transition of Si-buming into NSE 
is much smoother than the transition of C-O-burning into NSE. Furthermore, due to the more massive Si-S and C-O- 
layers for the Fe-core progenitors, the transitions take more time of the order 100 ms to fall onto the SAS. The presence 
of neutrino heating becomes important for the more massive Fe-core progenitors to drive the SAS to large radii on a 
longer timescale. The effects of the shock propagation across the transition between different thermonuclear regimes 
has been pointed out in Bruenn et al. (2006) with respect to the explosion dynamics in axially-symmetric simulations of 
massive Fe-core progenitors. In our spherically symmetric model we can not confirm the driving force of explosions of 
Fe-core progenitors to be the shock propagation across different thermonuclear regimes. It is rather the neutrino energy 
deposition, although enhanced, which initiates the explosion. The shock is accelerated additionally when crossing 
different thermonuclear regimes due to the density jumps at the transitions. 
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Figure 5: Heating (> 0) and cooling (< 0) rates for the 8.8 M progenitor model from Nomoto (1983,1984,1987) 
during the explosion phase at 20 ms (left panel), 25 ms (middle panel) and 30 ms (right panel) after bounce. For a 
better comparison, the velocities are plotted for the same post bounce times. 

3.3 Comparison of the neutrino spectra 

Striking is the agreement in the mean neutrino energies between all different progenitor models (including the O-Ne- 
Mg-core and the Fe-core progenitors) during the explosion phase, although the neutrino emissivities and opacities 
are enhanced for the Fe-core progenitor models (see Fig. [2]). The explosion phase for the O-Ne-Mg-core lasts only 
until about 40 ms post bounce and is significantly shorter in comparison to the more massive Fe-core progenitors. 
Furthermore, the luminosities are also smaller by a factor of 2. For all models, the electron antineutrino luminosity 
is larger than the electron neutrino luminosity on a timescale of 200 ms after the explosion has been launched. This 
slight difference reduces again at later times where the electron neutrino luminosity becomes again larger than the 
electron antineutrino luminosity. However, after the explosion has been launched the behavior of the luminosities are 
in qualitative agreement for all models. The same holds for the mean neutrino energies which increase continuously 
during the neutrino heating phase. The electron (anti)neutrinos have rms energies of about 12 (14) MeV where as after 
the explosion has been launched, rms energies of about 11 (13) MeV are obtained. The values remain constant on 
the timescale of 1 second post bounce. The (yu/r)-neutrinos have rms energies of about 18 MeV during the neutrino 
heating phase and about 15 MeV after the explosion has been launched. These differences in the mean neutrino energies 
and luminosities during the neutrino heating, initial and proceeding explosion phases are in correspondence with the 
electron fraction of the material, as will be illustrated in the following section. 

3.4 The electron fraction of the early ejecta 

During the neutrino heating phase, the neutrino spectra are mainly determined by mass accretion at the neutrinospheres. 
Neutron rich nuclei from the progenitor star with an electron fraction of Y e 0.45 are falling onto the oscillating SAS 
and dissociate into free nucleons and light nuclei, see Fig. [6] (d). These nucleons accrete than slowly onto the PNS 
surface at the center. Due to the increased electron-degeneracy behind the SAS in Fig. [6] (c), weak-equilibrium is 



established at a lower value of the electron fraction of Y e < 0.15. 
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Figure 6: Selected hydrodynamic variables during the initial explosion phase at three different post bounce times for 
the 10.8 M Q progenitor model from Woosley et al. (2002 ). 



As soon as the SAS is revived and propagates outward, see the velocity and density profiles in Fig. [6] (a) and (b), 
the electron degeneracy behind the expanding shock is reduced and weak-equilibrium is established at a larger value 
of the electron fraction of Y e > 0.56 [j The capture rates for electron-neutrinos at neutrons are favored over electron- 
antineutrino captures at protons. This slight difference results in an electron and hence proton excess. Consequently 
the explosion ejecta are found to be initially proton-rich. This behavior of the electron fraction was found for all our 
explosion models, for the 10.8 and 18 M Fe-core progenitors with artificially enhanced opacities and for the O-Ne- 
Mg-core using the standard opacities. Such explosion models were investigated with respect to the nucleosynthesis in 
general and with respect to the vp-process by Frohlich et al. (2006a-c). 

One of the main goals of the present investigation is to determine the behavior of the electron fraction for the 
initially proton-rich ejecta on a long timescale of the order 10 seconds, in a consistent manner. We explore the question 
if the material ejected in the neutrino driven wind becomes neutron-rich and which are the conditions (e.g. entropy per 
baryon, and expansion timescale) obtained in the neutrino driven wind. These aspects are of special relevance for the 
composition of the ejecta, which is determined via explosive nucleosynthesis analysis, in particular in order to be able 
to draw conclusions with respect to a possible site for the production of heavy elements via the r-process. Therefore, the 
continued expansion of the explosion ejecta must be simulated, for which the inclusion of a large physical domain of 
the progenitor up to the He-layer is required. Furthermore, since only the electronic charged current reactions including 
the neutrino fluxes can determine the electron fraction, the PNS contraction at the center and hence the contraction of 



3 The EoS from |Shen et al.|{l998| is limited to a maximum electron fraction of Y e > 0.564. The EoS has been extended by |Gogelein"|j2007) 
to model asymmetric nuclear matter with an electron fraction above 0.564. We assume an ideal gas of free nucleons and light nuclei, which is a 
sufficient assumption for the conditions found in the region of the extremely proton-rich ejecta. 



the neutrinospheres are essential. 



4 The neutrino driven wind 

In this section we investigate the post explosion evolution with special focus on the properties of the ejecta, in particular 
the electron fraction. We explore the problem if the initially proton-rich ejecta become neutron rich at later times of 
the order 10 seconds and if the conditions might indicate a possible site for the nucleosynthesis of heavy nuclei via the 
r-process. This has been assumed in static steady-state as well as parametrized dynamic wind models, based on the 
results obtained by Woosley et al. ( 1994} . The ejected material in their simulations does never become proton-rich, 



the electron fraction was found to continuously decrease with time after the explosion has been launched. Thus, the 
ejecta were investigated in a region where the conditions are favorable for the r-process. The stellar models applied 
to the present investigation of the neutrino driven wind are the 8.8 M Q O-Ne-Mg-core and the 10.8 and 18 M Fe- 
core progenitors, where for the latter two cases the explosions are obtained using the artificially enhanced opacities as 
described in §2.3. 

After the explosion has been launched, the region between the expanding explosion shock and the neutrinospheres 
cools rapidly and the density decreases continuously as illustrated in Figs. [7] and [8] (d) and (h). In order to determine 
the evolution of the electron fraction Y e , the non-local neutrino fluxes are as important as the local neutrino reaction 
rates Since the PNS and hence the neutrinospheres contract continuously due to the deleptonization, the degeneracy 
increases and matter at the PNS surface is found to be extremely neutron-rich with Y e < 0.1 (see Figs. [JJ and [HJ (d) and 
(f>. 

Independent of the progenitor model, the region on top of the PNS surface is continuously subject to neutrino 
heating during the post explosion phase. The dominant heating sources are the absorption of electron-(anti)neutrinos 
on free nucleons, due to the large fraction of free nucleons (dissociated nuclear matter) in the region on top of the 
PNS, as shown in Fig. [9] The neutrino pair production and thermalization processes have a negligible contribution 
to the heating outside the neutrinospheres. In order to compare the heating and cooling rates in Fig. [9] we plot the 
quantities with respect to the baryon density. While neutrino cooling is still dominantly present at ~ 500 ms post 
bounce (thin lines in Fig.|9]l, at later times after ~ 1 second post bounce (thick lines in Fig. |9|> neutrino cooling vanishes 
and only heating is found in the density domain of interest, i.e. between 10 7 - 10 12 g/cm 3 . Figs.|7](d) and[8](d) show the 
conditions for the contracting PNSs at the center via the radial baryon density profiles and the electron-neutrinospheres. 
The region of interest where the neutrino driven wind develops corresponds to the density domain of 10 7 - 10 10 g/cm 3 . 
The degeneracy of the early ejecta favors proton-rich matter where a large electron fraction of Y e ^ 0.54 is obtained. 
Hence, the absorption of electron-antineutrinos on free protons dominates over electron-neutrino absorption on free 
neutrons. The corresponding radial neutrino luminosities and rms energies are shown in Figs. [7] and [8] (e) and (g). 
In addition, for the first time we were able to follow the deleptonization burst from core bounce for several seconds 
over a large physical domain including several 10 s km of the progenitor star. The outrunning luminosity peak can be 
identified in the luminosities in Figs.[7]and[8](e) at 0.6 seconds after bounce at a distance between 5 x 10 4 and 10 5 km, 
leaving the computational domain between 1-2 seconds post bounce. 

After the explosion has been launched, the continued energy transfer from the neutrino radiation field into the fluid 
outside the neutrinospheres as illustrated in Fig. [9]drives the matter entropies to large values, shown in Figs. [7] and [8] 
(c). The heat deposition at the PNS surface accelerates matter to positive velocities, see Figs. [JJ and [8] (a), after ^ 1 
second post bounce. This matter outflow is known as the neutrino driven wind, which proceeds adiabatically at larger 
radii. This is consistent with the constant radial entropy per baryon profiles in the graphs (c). Furthermore, the rapidly 
decreasing luminosities reach values below 5 x 10 51 erg/s already 1 second after bounce (see Figs. [7] and [8] (e)). The 
luminosities continue to decrease and reach values below 1 x 10 51 erg/s at 10 seconds after bounce. The mean neutrino 
energies also decrease constantly where values below 10 MeV for the electron-flavor neutrinos and below 12 MeV for 
the (/z/r)-neutrinos are obtained (see Figs. [7J and [H](g)). 

Several previous wind studies achieved supersonic matter outflow velocities for the neutrino driven wind due to 
assumed large luminosities. In any case, the accelerated material of the neutrino driven wind collides with the slower 
and subsonically expanding explosion ejecta. Comparing Figs. |7j and [8] the more compact wind region of the 18 M Q 
progenitor model produces a less pronounced neutrino driven wind in comparison to the 10.8 M progenitor model. 
The densities of the wind region are larger up to two orders of magnitude and the temperatures are higher by a factor 
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Figure 7: Selected hydrodynamic variables during the formation of the neutrino driven wind at three different post 
bounce times for the 10.8 M progenitor model from Woosley et al. (2002 1. In addition, graphs (e) and (g) show the 
neutrino luminosities and rms neutrino energies (solid lines: v e , dashed lines: v e , dash-dotted lines: v^/ T ). For this 
progenitor model the neutrino driven wind becomes supersonic, using the enhanced opacities. 
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Figure 8: The same configuration as Fig.[7]for the 18 M progenitor model from Woosley et al. ( 2002) 1. The neutrino 
driven wind remains subsonic for this progenitor model, even using the enhanced opacities. 
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Figure 9: Heating and cooling rates with respect to the baryon density at three different post bounce times 0.5 second 
(thin lines), 1 second (intermediate lines), 5 seconds (thick lines) for the 10.8 M progenitor model from Woosley et al. 
( |2002| . 



of 2. The resulting velocities of the neutrino driven wind outflow are smaller by a factor of 2 and stay below 10 4 
km/s. Hence, the neutrino driven wind remains subsonic for all times for the 18 M Q progenitor model (see Fig. [8] 
(a)) where the wind develops supersonic velocities for the 10.8 M Q progenitor model (see Fig. [7] (a)). In the case 
of a supersonic neutrino driven wind, this leads to the formation of a reverse shock known as the wind termination 
shock. The formation of the wind termination shock of the 10.8 M progenitor model is illustrated in Fig. 10 and 
will be discussed in the following paragraph. In order to analyze the dynamical evolution and the consequences of the 
formation of the reverse shock, steady-state wind models can not be used. Radiation hydrodynamics is required in order 
to describe the dynamical effects consistently. Our results, obtained using general relativistic radiation hydrodynamics 
based on spectral three flavor Boltzmann neutrino transport, are in qualitative agreement with the detailed parametrized 
investigation by Arcones et al.| ( [2007] ). 



During the initial and subsonic wind expansion, the matter entropies in Fig. 10 (c) increase slowly from 4 to 5 - 10 
kB/baryon and the densities in Fig. 10 (b) and temperatures in Fig. 10 (e) decrease on a long timescale over several 
seconds. Furthermore, the reduced degeneracy in the wind increases the electron fraction shown in Fig.[l0](d) slowly 
on the same timescale. When the material is accelerated supersonically with velocities of several 10 4 km/s up to radii of 
a several 10 3 km (see Figs. 



10 



' (a) and (f)), the entropies increase from s ^ 5-10 kB/baryon to s ^ 40 - 60 kB/baryon on 
a short timescale of the order of 100 ms. During this rapid expansion, the density and temperature decrease drastically 
from 10 10 g/cm 3 to 10 4 - 10 2 g/cm 3 and from 3 MeV to 0.1 - 0.01 MeV respectively (see Figs.fTol(b) and (e)). It also 



corresponds to a rapid decrease of the degeneracy which in turn is reflected in a rapid increase of the electron fraction 
of the accelerated material on top of the PNS surface, from Y e 0.1 to Y e 0.56 (see Fig. 10 (d)). Furthermore, the 
supersonically expanding neutrino driven wind collides with the explosion ejecta as can bee seen in Fig. 10 (a) (solid 
red line) at radii of several 10 4 km. Consequently, the previously accelerated material decelerates behind the explosion 
ejecta as can be seen in the velocities in Fig. 10 (f). This phenomenon becomes significant after about 2 seconds 
post bounce and corresponds to the formation of the reverse shock, i.e. the wind termination shock, (see Fig. [T0|(a) 
dashed red line at radii of several 10 3 km). It causes an additional entropy increase to the final values of s ^ 50 - 100 
k B /baryon. During the rapid deceleration on the same short timescale of the order 100 ms, the densities in Fig. 10 (b) 
and temperatures in Fig. 10 (e) increase again slightly, where the degeneracy increases and hence the electron fraction 
reduces slightly to values of Y e ^ 0.54. The following dynamical evolution is given by the subsonic and adiabatic 
expansion of the explosion ejecta on a longer timescale of the order of seconds. The density and temperature decrease 
slowly where the entropies of s ^ 50 - 100 kB/baryon and the electron fraction of about Y e = 0.54 remain constant. 
The latter aspects are essential for the nucleosynthesis analysis of the ejecta. It can be understood in the sense that the 
neutrino reaction rates freeze out and the matter conditions correspond to the neutrino free streaming regime. 

Note that the strong neutrino driven wind for the 10.8 M progenitor model is obtained using the enhanced opacities 
as introduced in §2.3. We additionally illustrate selected properties of the neutrino driven wind for the 8.8 M Q progen- 



itor model in Fig. 11 where a strong neutrino driven wind was obtained using the standard emissivities and opacities 



given in |Bruenn| ( |1985| l. This is due to the low density of the region between the neutrinospheres at the PNS surface 
and the expanding explosion shock, where neutrino heating via the standard rates and energy from nuclear burning are 
sufficient to drive a strong supersonic matter outflow. Matter entropies increase to s ^ 10 kfi/baryon during the initial 
acceleration of the wind and the densities and temperatures decrease slowly on a timescale of seconds. The properties 
during the initial acceleration observed are similar to those of the more massive 10.8 M Q Fe-core progenitor. The same 
holds for the acceleration to supersonic velocities. The timescale is reduced to 100 ms where the entropies increase 



rapidly to s ^ 20 - 50 ke/baryon (see Fig. 1 1 (c)) and due to the reduced degeneracy the electron fraction increases 
from Y e = 0.1 at the PNS surface to Y e = 0.56 (see Fig. 11 (d)). Density and temperature decrease to 10 - 100 g/cm 3 



and 0.001 MeV respectively (see Fig.[TT](b) and (e)). The difference to the more massive 10.8 M Q Fe-core progenitor 
is due to the smaller mass enclosed between the PNS surface and the expanding explosion ejecta. For the more massive 
10.8 M Fe-core progenitor in Fig. 10 (f), the previously accelerated material collides with the explosion ejecta already 
after a few 100 ms. Here the supersonic wind expands on a much longer timescale up to several seconds before it col- 
lides with the explosion ejecta (see Fig. 1 1 (f)). During this adiabatic expansion, entropy and electron fraction remain 
constant. The fast material collides with the much slower expanding explosion ejecta so that the material is decelerated 
and the reverse shock appears. This is again similar to the formation of the reverse shock for the more massive 10.8 
M Q Fe-core progenitor as discussed above. Matter entropies increase to several 100 kB/baryon (see Fig{TT|(c)), density 
and temperature increase sightly (see FigJTTl(b) and (e)) and the electron fraction reduces slightly to Y e ^ 0.52 - 0.54 
due to the increased degeneracy (see Figjl l|(d)). The following evolution is determined by the adiabatic expansion of 
the explosion ejecta during which the entropy and electron fraction remain constant. 

In the following paragraph, we will discuss the composition of the neutrino driven wind region to some extent. This 
is possible due to the recently implemented nuclear reaction network. It includes the free nucleons and the symmetric 
nuclei from 4 He to 56 Ni plus 53 Fe, 54 Fe and 56 Fe. The initial composition is given by the progenitor model. Mostly 
28 Si and 3l) S are shock-heated and burned to Fe-group nuclei due to the temperature and density jump during the initial 
expansion of the explosion shock (see Fig. [12] and compare with Figs. [JJ and [8] (d) and (f)). The large fraction of 
these Fe-group nuclei reduces behind the explosion shock due to photodisintegration, indicated by the region of low 
density and high entropy in Figs. [7J and [8] (b) and (c). This produces a large fraction of a-particles, which in our 
model represent light nuclei. The region of a-particle domination behind the expanding explosion shock increases 
with time. This behaviour is illustrated in Fig. 12 for both Fe-core progenitor models under investigation. The position 
of the explosion shock coincides with the maximum of the mass fraction of Fe-group nuclei (in particular 56 Ni). In 
addition, density and temperature of the neutrino driven wind on top of the PNS surface decrease continuously with 
time. The low temperatures and densities in that region do not justify the assumption of NSE beyond ~ 1 second after 
bounce, where temperatures reach values below 0.5 MeV. Instead, our nuclear reaction network is used to determine 
the composition in that region. The decreasing density and temperature and the presence of a large fraction of free 
nucleons favor the freeze out of light nuclei. Finally, the entire region between the expanding explosion shock and the 



PNS surface is found to be dominated in our simulations by <?-particles. In Fig. 12 the radii of the NSE to non-NSE 
transitions are indicated by vertical lines. The slight mismatch between the abundances between the heavy 'Fe'-group 
nuclei (the representative heavy nucleus with average atomic mass and charge in NSE) and 56 Ni (non-NSE) as well as 
between the a's is due to the different nuclear models used for the two regimes. While in NSE the EoS for hot and dense 
nuclear matter assumes 56 Fe as the most stable nucleus due to the smallest mass per nucleon for low temperatures and 
densities, the nuclear reaction network applied in non-NSE calculates the composition dynamically based on tabulated 
reaction rates. 



5 Comparison with previous wind studies 



5.1 The proton-to-baryon ratio of the wind 

The approximations made in previous wind studies is a simplified description of the radiation-hydrodynamics equa- 
tions, see for example Duncan et al. (1986) and Qia n and Woosley] ( |1996 1. More crucial is the absence of neutrino 
transport. Neutrino heating and cooling is calculated only based on parametrized neutrino luminosities and mean en- 
ergies. Hence, such models explore the neutrino driven wind by varying the neutrino luminosities and energies, where 




Figure 10: Evolution of selected mass shells in the neutrino driven wind (as listed in graph (a) from 1.4429 - 1.4437 
M Q baryon mass) for the 10.8 Mq progenitor model from Woosley et al. (2002 ) where the enhanced opacities are used. 
Graph (a) shows in addition gray-scaled the entropy per baryon, the position of the expanding explosion shock (thick 
solid line) and the position of the wind termination shock (thick dashed line). 
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Figure 11: Evolution of selected mass shells in the neutrino driven wind (as listed in graph (a) from 1.34718 - 1.34750 
M Q baryon mass) for the 8.8 M progenitor model from Nomoto (1983,1984,1987) where the standard emissivities 
and opacities given in Bruenn ( 1985 ) are used. The graphs show the same configurations as Fig.fTO 
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Figure 12: Selected radial mass fraction profiles during the neutrino driven wind phase for the 10.8 M Q (left panel) 
and the 18 M (right panel) progenitor model from Woosley et al. (2002| at 1 second (top), 3 seconds (middle) and 
5 seconds (bottom) post bounce. The vertical lines represent the separation of NSE (EoS for hot and dense nuclear 
matter) where heavy nuclei are represented by a single Fe-group nucleus 'Fe' with average atomic mass and charge 
and non-NSE (nuclear reaction network) where the most abundant Fe-group element is 56 Ni, at temperatures of ^ 0.5 
MeV. 



the simplified radiation-hydrodynamics equations are solved (see for example Thompson et al. ( 2001| l). Since neutrino 
transport is neglected, the evolution equation for the electron fraction Eq. (|15j> cannot be solved consistently because 
the neutrino distribution functions are unknown. In the following paragraph, we will discuss the assumptions made for 
the evolution of the electron fraction in the neutrino driven wind which go back to Qian and Woosley ( 1996 ). 

Applying the theory of weak interactions based on the reaction rates A,j for the reaction particles (z, J), i.e. electron 
and positron as well as electron neutrino and antineutrino captures, the evolution equation for the electron and positron 
fractions can be written as follows 

dY.- 

A Ven Y v Y n , (16) 



dY e - 
~~dT 
dY e * 
~dT 



(17) 



These expressions can be combined to calculate the evolution of the total number of charges, using the relations Y p 
and Y n = 1 — Y e , 



dY e 
dt 



= - (Y e - - Y e+ ) 
dt 



\A e - p Y e * + A e + n Y e + + A Ve „Y Ve + A Ve pYy e ^ Yg. 



assuming fully dissociated nuclear matter. This expression is approximated in a crucial but powerful way (IQian and 



Woosley (1996) Eq.(73)), ignoring contributions from electron and positron captures as well as the decoupling of 



radiation and the angular dependency of the neutrino distribution function on the distance from the energy-dependent 
neutrinospheres, as follows 



Y„ * 



A-e-pYe* + A e +„Y e + + A Ve „Y Vi 



+ A-V e pYv, 



A-v e nYv e / ^v e pYy e 



(18) 
(19) 



This approximation of the electron fraction was further simplified and expressed in terms of the neutrino luminosities 
L v and mean neutrino energies (E v ) and the well known rest mass difference between neutron and proton Q - m n -m p = 
1.2935 MeV, as follows 

i.2 q 2 y 1 



Y„ a 



1 + 



, (E 9 ,)-2Q + 



(20) 



which is used in previous static steady-state and parametrized dynamic studies of the neutrino driven wind. 
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(a) Different electron fraction approximations (dash-dotted 
line: Eq. |19| dashed line: Eq. |20) in comparison with Boltz- 
mann transport (solid line). 
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(b) Evolution of the electron fraction according to Eq. |20j. 
The three lines correspond to different assumptions for the 
electron antineutrino mean energy. 



Figure 13 

10.8 M progenitor model from Woosley et al 



The Electron fraction approximations at a distance of 10 km outside the electron-neutrinosphere for the 

2002) . 



Fig. 13(a) compares the electron fraction behavior at a distance of 10 km outside the electron-neutrinosphere, 



from Boltzmann neutrino transport (blue solid line) with the approximations Eq. ( 19 1 based on the neutrino capture 
rates (red dashed line) and Eq. ( f2"0) i based on the luminosities and mean neutrino energies (red dash-dotted lines). The 
approximations are in qualitative agreement with Boltzmann transport. The differences on the longer timescale are most 
likely due to the presence of light and heavy nuclei which are not taken into account explicitly in the approximations. 

All descriptions agree 



17) 



They change the number of free nucleons available for the reactions in Eq. (16i and 
qualitatively in the prediction of a generally proton-rich material in the wind, based on the neutrino spectra obtained 
via Boltzmann transport. 



5.2 The neutrino observables in the wind 



Comparing the neutrino spectra in Fig. 14 with the spectra assumed in previous static steady-state and dynamic wind 
studies (see for example Thompson et al. (2001 ) and Arcones et al. (2007i), we find two major differences: One, 
the neutrino luminosities and mean neutrino energies assumed are significantly larger than those we find and two, the 
assumed behaviour with respect to time is different. 

The commonly used assumptions made in static steady-state and parametrized dynamic wind studies go back to the 
detailed investigation from Woosley et al. ( 1994[ ), who performed core collapse simulations based on sophisticated input 



physics. They investigated the neutrino driven explosion of a 20 M progenitor star and followed the evolution for 18 
seconds post bounce into the neutrino driven wind phase. In their simulations the electron (anti)neutrino luminosities 
decreased from initially 4 x 10 52 (3 x 10 52 ) erg/s at the onset of the explosion to 6 x 10 51 (7.5 x 10 50 ) erg/s at 10 
seconds after bounce, where strictly Ly e > L Ve after the onset of the explosion. The difference between the neutrino and 
antineutrino luminosities remained small and constant with respect to time up to 3 seconds post bounce and increased 
only significantly after 4-5 seconds post bounce, after which the difference reached its maximum of 1.5 x 10 50 
erg/s at the end of the simulation at about 18 seconds post bounce. The electron flavor neutrino luminosities in our 
models follow a different behaviour. They reach 1 x 10 51 erg/s at about 5, 6 and 8 seconds post bounce for the 8.8, 
10.8 and 18 M progenitor models respectively. The larger electron flavor neutrino luminosities for the more massive 
progenitors are in correlation with the more massive PNSs and the hence larger number of neutrinos emitted. However, 
the difference between electron-neutrino and electron-antineutrino luminosities found in the present investigation is 



significantly smaller than the difference in Woosley et al. ( 19941. During the initial explosion phase until about 300 



ms after the onset of the explosion, the electron antineutrino luminosity is slightly larger than the electron neutrino 
luminosity by about 1 x 10 50 erg/s which explains the electron fraction of Y e > 0.5 of the early explosion ejecta. After 
about 900 ms post bounce, the luminosities can hardly be distinguished where during the initial neutrino driven wind 
phase after about 1 second after bounce the electron neutrino luminosity becomes larger than the electron antineutrino 
luminosity by about 1 x 10 50 erg/s. This difference reduces again at later times at about 6 seconds post bounce and the 



electron flavor neutrino luminosities become more and more similar (see Fig. 14 1 



Even more different are the values and the behavior of the mean neutrino energies, see Fig. 14 and compare with 
Fig. 2 of Woosley et al. ( 1994| >. They found (p/r)-neutrino energies of about 35 MeV which remained constant with 
respect to time. Their electron-antineutrino energies increased slightly from about 20 MeV to 22 MeV where the 
electron-neutrino energies decrease from 14 MeV to 12 MeV. This increasing difference between the electron neutrino 
and antineutrino spectra favored neutron-rich material, which was consistent with their findings of Y e < 0.5 for the 
material ejected in the neutrino driven wind in Woosley et al. ( 1994) >. We cannot confirm these results for the mean 
neutrino energies nor the evolution of the spectra. In contrast, all mean neutrino energies decrease with respect to time 
for all our models. This is a consequence of lepton number and energy loss of the central PNS where the neutrinos 
diffuse out. The electron (anti)neutrino energies decrease from about 10 (12) MeV at the onset of the explosion to 
about 8.5 (9) MeV and the (yu/T)-neutrino energies decrease from 16 MeV to 10 MeV at the end of the simulations. 
Hence, not only the mean energies decrease also the difference between the electron flavor neutrino spectra decreases. 
The reason for the neutrino spectra to become more similar with respect to time is related to the evolution of the 
thermodynamic properties at the neutrinospheres, and will be discussed in the following subsection. 



5.3 The PNS contraction 

The behaviour of the neutrino spectra and hence the evolution and the properties of the neutrinospheres is related to the 
PNS contraction. The contraction is caused by a continuous deleptonization and translates to a continued steepening of 
the density gradient at the PNS surface. Hence, the neutrinosphere radii for the electron flavor neutrinos move closer 
together with time. The evolution of the neutrinosphere radii for both electron neutrino and antineutrino are illustrated 
in Fig.[l5](a) for the 10.8 M Q progenitor model. Their difference reduces from 740 m at about 1 second post bounce to 
370 m at about 5 seconds post bounce and further to 260 m at about 10 seconds post bounce. 

This contraction behaviour has consequences for the neutrino spectra, which are determined during the neutrino 
driven wind phase by diffusion rather than by mass accretion. Hence, the electron flavor neutrino luminosities can be 
determined as follows 



1 



-4nr z 



(21) 



where u Y oc T is the thermal black body spectrum for ultra-relativistic fermions with temperature T. The matter 



temperatures at the neutrinospheres decrease with respect to time as shown in Fig. 15 (b). This is due to the continued 



loss of lepton number and energy, carried away by the diffusing neutrino radiation field as illustrated in Fig. 16 for 
the 10.8 M Q from Woosley et al. ( 2002| l. The lepton number decreases from Yi - 0.3 at 2 seconds post bounce to 
Yi =i 0.18 at 10 seconds post bounce, see Fig. 16 (c). The additionally reduced mean neutrino energies (on average), 
from {E) lms 150 MeV to (E) lms 50 MeV (see Fig. [16] (d)), and the consequent reduced temperature-peak inside 



the PNS (see Fig. 16 (b)) cause the contraction of the outer layers of the PNS. This can be identified via the density 
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Figure 14: Neutrino luminosities and mean energies with respect to time after bounce for the 8.8 M Q O-Ne-Mg- 
core from Nomoto (1983,1984,1987) (left panel) and the 10.8 M (middle panels) and 18 Mq (right panel) Fe-core 



progenitor models from Woosley et al. ( 2002| l, measured in a co-moving reference frame at a distance of 500 km. 



increase in Fig. 16 (a). Note in addition to the maximum temperature decrease at the outer layers of the PNS, the 
central temperature increases from 18 MeV to 23 MeV on the post bounce times between 2 and 10 seconds. This 
is caused by the contraction of the deleptonizing outer layers of the PNS which compresses the central part. The 



evolution of the reducing temperature at the neutrinospheres is shown in Fig. 15 (b). In combination with the loss of 
leptons number, it explains the decreasing electron flavor neutrino luminosities and mean neutrino energies with respect 
to time. Furthermore, the temperature difference decreases with respect to time from 0.467 MeV at about 1 second 
post bounce to 0.362 MeV at about 10 seconds post bounce. Consequently the neutrino spectra become more similar 
with respect to time. As illustrated in Fig. 14 the difference in the electron flavor neutrino luminosities and mean 



neutrino energies decreases for all models under investigation. In addition, Fig. 15 (c) and (d) illustrate the evolution of 
the baryon density and the electron fraction at the corresponding neutrinospheres. It becomes additionally clear from 



the electron fraction approximation Eq.( 20 1, that it is not the absolute values for the mean neutrino and antineutrino 



energies that determine whether matter becomes neutron- or proton-rich but their ratio. 
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Figure 15: Evolution of the neutrinosphere radii in graph (a) and temperature and density at the corresponding neutri- 
nospheres in graphs (b) and (c) respectively for the 10.8 M Q progenitor model from Woosley et al. (2002 ). 
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Figure 16: Radial PNS profiles at two different post bounce times (thick lines: 2 seconds, thin lines: 10 seconds), for 
the 10.8 M Q progenitor model from Woosle y"et al.| ( |2002| . 



Since this difference is small in our simulations, with initially (E Ve ) rms ^10 MeV and {EyJ}^ ^13 MeV (at about 
1 second post bounce) and only (E Vt ) - 9 MeV and (E 9e ) rms ^11 MeV (at later times at 10 seconds post bounce), 
the values found for the electron fraction of Y e > 0.5 (solid line in Fig. 13(b) for the ^-approximation based on the 
luminosity and mean neutrino energies) clearly illustrate that the accelerated matter in the neutrino driven wind stays 
proton-rich for more than 10 seconds. This is in qualitative agreement with Boltzmann transport as discussed above 
and shown in Fig. |13(a)| Hence we find Eq.( [20] > to be a good approximation to model the electron fraction in the 
wind. On the other hand, most of the previous studies select the neutrino luminosities and mean energies to investigate 
a neutron-rich neutrino driven wind. In order to test the appearance of Y e < 0.5 under such conditions, we increase 
the difference between the mean neutrino and antineutrino energies by hand. We evaluate expression ( 20 1 shown in 



Fig. 13(b) at 10 km outside the electron-neutrinosphere for 1.2 (dashed line) and 1.5 (dash-dotted line) times larger 



electron-antineutrino mean energies. For the first value, Y e decreases but matter remains slightly proton-rich, which 
for the latter value matter becomes neutron-rich. Indeed, the larger the difference between neutrino and antineutrino 
spectra are, the smaller becomes the electron fraction in the wind. Note that the luminosities and electron-neutrino 
energies remained unmodified for this experiment. Such an increase of the energy difference between neutrinos and 
antineutrinos could perhaps be related to the uncertainty of the EoS for nuclear matter, which will be discussed in the 
following paragraph. 

The assumed PNS radii in previous wind studies reach about 10 km shortly (< 1 second) after the explosion has 
been launched. We define the radius of the PNS as the position of the electron-neutrinosphere at the steep density 
gradient at the PNS surface. The approximated inner boundary of the physical domain in most wind models is close 
to but still inside this radius. The position of the neutrinospheres and the contraction of the PNSs found in the present 
paper differ significantly from the assumptions made in most previous wind studies. We find PNS radii of about 40 km 
at the time of the explosion and 20 km at about 2 seconds after bounce. During the later evolution, the PNS contraction 
slows down. The PNS profile and hence the position of the neutrinospheres as well as the contraction behavior itself 
is given implicitly by the EoS for hot and dense nuclear matter as well as the PNS deleptonization. For the stiff EoS 
from |Shen et al. ( 1998] ) and both the 10.8 and 18 M Q progenitors, the PNSs reach radii of 14.5 - 15 km only at about 
10 seconds after bounce (see Fig. 15 1. The larger radii of the neutrinospheres result in lower neutrino luminosities and 
mean energies and a smaller difference between neutrino and antineutrino spectra in comparison to the assumptions 
made in most previous wind models. This is in agreement with Arcones et al. ( 2007[ > who additionally assume PNS 
radii of 15 km and find conditions that differ more from previous wind studies. They obtained significantly larger 
values for the electron fraction. To summarize, this effects and the different behaviour of the neutrino spectra assumed 
in the previous wind studies leads to different matter properties of the neutrino driven wind. A detailed comparison 
study of fast and slowly contracting PNSs with respect to the neutrino driven wind, e.g. applying EoSs with different 
compressibilities and asymmetry energies, would be necessary in the context of radiation hydrodynamics simulations 
using spectral Boltzmann neutrino transport. 



6 Long term post bounce evolution 

During the neutrino driven wind phase, the neutrino luminosities and mean neutrino energies decrease continuously, 



which leads to a constant decrease in the net-heating rates. At luminosities below 10 51 erg/s (see Fig. 14 1, the supersonic 
matter outflow for the 10.8 M progenitor model descends into a subsonic expansion. The wind termination shock turns 
again into a subsonic neutrino driven wind. At later times, the neutrino driven wind settles down to a quasi-stationary 



state with no significant matter outflow, illustrated at the example of the 18 M Q progenitor model in Fig 17 (a). The 
explosion shock continues to expand and the material enclosed inside the mass cut accretes onto the PNS at the center. 
In combination with the deleptonization, this leads to the continuous PNS contraction. However, the contraction 
proceeds on a timescale of seconds and hence the PNS can be considered in a quasi-stationary state. The dense and 
still hot and lepton-rich PNS at the center is surrounded by a low density and high entropy atmosphere, composed of 
light and heavy nuclei. See for example the abundances of the 18 M progenitor for the post bounce time of 22 seconds 



in Fig. 17 (e). 

The internal temperature profile of the PNS is not constant. The central region of the PNS did not experience 
shock heating immediately after the Fe-core bounce, since the initial shock forms at the edge of the bouncing core. 
Its mass scales roughly with Y% and is typically around values of 0.5 - 0.6 M for low and intermediate mass Fe-core 
progenitors. Hence, the central temperature after bounce is given by the thermodynamic conditions at bounce. The 
temperature changes only during the post bounce evolution due to compressional heating and the diffusion of neutrinos. 
The shock heated material inside the PNS shows significantly higher temperatures than at the center. The temperature 
decreases again towards the PNS surface where the matter is less dense (for the illustration of the radial temperature 



profile inside the PNS as well as the dynamical evolution of temperature and density, see Fig 17 (f) and (b) at selected 
post bounce times between 5-22 seconds). The neutrinos diffuse continuously out of the PNS and carry away energy. 
The central electron fraction reduces from Y e 0.25 at the onset of the explosion to Y e 0.15 at 22 seconds after 



bounce (see Fig 17 (d)). It relates to a temperature decrease from about 35 MeV initially (at 3 seconds post bounce) 



to 23 MeV at about 22 seconds post bounce. This corresponds to the initial and neutrino dominated cooling phase. 
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Figure 17: Radial profiles of se- 
lected hydrodynamic variables 
for the 18 M Q progenitor model 
at three different post bounce 
times, illustrating the disappear- 
ance of the neutrino driven wind 
and the PNS cooling and con- 
traction. Graph (e) illustrates the 
composition at 22 seconds post 



Radius, r [km] bounce. 



Unfortunately the achieved temperatures are not representative since important neutrino reactions, such as the direct 
and modified Urea processes, are not yet taken into account. 




Figure 18: Comparing selected hydrodynamic variables using the standard reactions rates based on |Bruenn ( 1985] ) 
(solid lines) and the artificially enhanced rates (dashed lines) for the 10.8 M (top) and the 18 M progenitor model 
(bottom). 



7 Discussion 



The neutrino driven wind was found to occur in all three progenitor models under investigation, the 8.8 M Q O-Ne- 
Mg-core and the 10.8 and 18 M Fe-core progenitor models. Because the neutrino driven explosions for the Fe-core 
progenitors are launched using artificially enhanced neutrino reaction rates, one may ask about the impact of these 
modified rates on the neutrino driven wind. Therefore we performed additional runs for which we switch back to 
the standard opacities given in Bruenn ( 1985] > after the explosion has been launched. The times when we switch 
back is about 500 ms after bounce, chosen such that the dynamics of the explosion ejecta does not change anymore 
significantly due to neutrino heating. However, the lower opacities translate to a significantly smaller net-heating by 
a factor of 5 - 6 in the region on top of the PNS where the neutrino driven wind develops. The energy deposition is 
still sufficient to drive the neutrino driven wind but the matter velocities are smaller by a factor of 2 - 5 in comparison 



to the wind velocities using the enhanced reaction rates (see Fig. 18 (a)). The main effect of the artificially enhanced 
reaction rates and the hence increased neutrino heating on the dynamics is clearly the stronger neutrino driven wind. 
For the 10.8 M progenitor model and with the enhanced heating, the wind even develops supersonic velocities (as 



discussed above in § 4) in Fig. 18 (a) (top panel). The supersonic wind collides with the explosion ejecta where matter 



decelerates and hence the reverse shock forms, which additionally increases the entropy in the wind (see Fig. 18 c) 



(top panel)). This additional entropy increase is absent in the simulations using the standard opacities, where the 



wind stays subsonic. This is also the case for the 18 progenitor model (Fig. 18 bottom panel), with and without 
the enhanced opacities. The neutrino driven wind of the O-Ne-Mg-core is illustrated in Fig. [TT] using the standard 
rates based on Bruenn (1985). The formation of a supersonic neutrino driven wind could be confirmed including the 
formation of the wind termination shock. Hence, one may speculate whether only low mass progenitors develop a 
strong neutrino driven outflow where the influence of the wind to the matter properties of the ejecta becomes small 



for more massive progenitors. The progenitor dependency of the neutrino driven wind is related to the density of the 
envelope surrounding the PNS after the explosion has been launched, which is significantly higher for more massive 
progenitors and hence the neutrino driven wind is weaker. 

The agreement of the time evolution of the mean neutrino energies between all three progenitor models under 
investigation (using the enhanced and standard opacities) in Fig. [14] is striking. The impact of the artificial heating 
to the neutrino observables and hence to the electron fraction in the wind is less pronounced. The influence on the 



composition of the wind is illustrated via the electron fraction in Fig. 18 (d). Using the standard rates, the wind stays 
slightly more proton-rich. Increasing the charged current reaction rates allows /3-equilibrium to be established on a 
shorter timescale. In addition, matter stays slightly more proton-rich for the the less intense neutrino driven wind, 
which develops for the Fe-core progenitors using the standard neutrino opacities. The additional electron fraction 
reduction in the neutrino driven wind for the models using the enhanced neutrino reactions is found due to the larger 
degeneracy obtained in the stronger deceleration behind the explosion ejecta. However, the findings of generally 
proton-rich ejecta as well as the generally proton-rich neutrino driven wind does not change. The corresponding 



densities and entropies per baryon in the wind are shown in Fig. 18 (b) and (c). The effects of the artificial heating are 



slightly smaller entropies. The higher matter outflow velocities in the wind region using the artificial heating results 



additionally in lower densities, shown in Fig. 18 (b). 

The artificially increased charged current reaction rates cannot be justified by physical uncertainties of the rates 



themselves. Similar to the large luminosities assumed in Arcones et al. (2007 1, they could rather be seen as a lowest 
order attempt to take the effects of multi-dimensional phenomena into account. For example, known fluid instabilities 
increase the neutrino energy deposition efficiency. Convection allows matter to stay for a longer time in the neutrino 
heating region (see Hera nt et al. | ( |1994| , |Janka and Muller~| ( |1996| )). Present axially symmetric core collapse models of 
massive Fe-core progenitor stars (even non-rotating) predict bipolar explosions (see Janka et al. ( 2008 )). The deviation 
from a spherical description and hence the deformation of the SAS due to fluid instabilities takes place during the 
neutrino heating phase on a timescale of several 100 ms after bounce. In multi-dimensional models, the luminosities 
are powered by a significantly larger mass accretion because the up-streaming neutrino heated matter is accompanied 
by a large down-stream of cold material. These larger luminosities may power a supersonic neutrino driven wind 
behind the explosion ejecta while the neutrino driven wind may remain absent in the angular directions of the accreting 
material which will not be ejected. 

Previous wind models have long been investigated as a possible site for the production of heavy elements via the 
r-process, motivated by the expectation of the ejection of neutron rich material, the large entropies in the neutrino 



driven wind and the short timescale of the neutrino driven wind expansion (see Hoffman et al. (2007) and Panov and 
Janka (2009} and references therein). The relevant quantities are shown in Figs. 10 and 11 for the 10.8 M Q and 8.8 
M Q progenitor models respectively. Illustrated are several selected mass shells that are part of the region where the 
neutrino driven wind develops in our radiation hydrodynamics model based on spectral three flavor Boltzmann neutrino 
transport. The inclusion of neutrino transport in a dynamical model is essential in order to obtain consistent neutrino 
spectra which determine the evolution of the electron fraction and the PNS contraction due to deleptonization and 
mass accretion. In comparison to previous static steady-state and dynamic wind models - where these ingredients 
were assumed - we confirm several properties of the accelerated material in the neutrino driven wind, such as the fast 
expansion and the large matter outflow rate shown in Fig. 19 a) and (b), the high velocities in the Figs. 10 f) and[TTJf)) 
and the rapidly decreasing density and temperature of the accelerated material in Fig. 10 (b) and (c) respectively. The 



expansion timescale in Fig. 19 b) is given by the following expression 



^"dyn 



r=().5MeV 



evaluated at the surface of constant temperature of T = 0.5 MeV, compared with an alternative definition of the 
dynamic timescale which has been introduced in Thompson et al. ( 2001] > 

- 1 



1 1 dp 
v p dr 



r=0.5MeV 



as well as a timescale approximation which has been derived in Qian and Woosley ( 1996 1 Eq.(61) 

1 1 



-/?p N sMpNS, 



where additionally an approximation for the mass outflow rate is derived as follows 



dM 
~dt 



QW 



K ,5/3 10/3 ,,5/3 j.-2 



where /?pns and Mpns are the PNS radius and mass respectively, which we take to be given by the electron-antineutrinosphere, 
and €y t = (e^ ^ / (Ey s ) where (E Ve ) is the mean electron-antineutrino energy and L Ve is the electron-antineutrino lumi- 
nosity, both taken at the surface of constant temperature of T — 0.5 MeV. In comparison with previous wind studies 

30- 100 ms for the 
(a) and (b) 
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(see for example Fig. 4 of lArcones et al. (2007 1), we find generally a longer timescale of Td yn 
neutrino driven wind which corresponds to a mass outflow rate of 10~ 3 - 10~ 5 M Q s _1 shown in Fig. 
respectively. However, we find generally a much longer timescale than the approximation tq W which is of the order 
~ 1 — 5 ms. The same holds for the mass outflow rate dMldt\qyf which is smaller by one order of magnitude. This may 
be related to the crucial assumptions made during the derivation of the above expressions, e.g. hydrostatic equilibrium, 
^pns = 10 km, L Ve =t L Ve , £y e = 20 MeV, which differ significantly from our findings. However, the wind entropies of 
40-100 kfi found (initially driven due to neutrino heating and additionally due to the deceleration in the reverse shock) 
are significantly smaller than often assumed in the literature and the previously accelerated matter does not become 
neutron rich as the neutrino driven wind decelerates behind the explosion ejecta but stays slightly proton-rich where 
Y e =i 0.54 for more than 10 seconds. This, in combination with the much slower PNS contraction illustrated via the 



neutrinospheres in Fig. 15 in comparison to static steady-state and dynamic wind models suggest that the assumptions 
made in previous wind studies should be carefully reconsidered. With respect to Woo sley et al.| ( [T994] l (e.g. Fig. 3), 
we find generally smaller mean neutrino energies which decrease with respect to time after bounce. This results in a 



decreasing difference between the electron flavor neutrino mean energies, while in Woosley et al. ( 1994) this difference 
increases. This fact in combination with the different PNS properties found in Woosley et al. ( 1994 1, enabled a strong 
neutrino driven wind where high entropies up to 400 kfi/baryon and a low electron fraction of Y e 0.35 - 0.45 was 
obtained. These properties of the neutrino driven wind differ quantitatively from our results, where only entropies of 
60 - 100 ke/baryon are obtained and matter stays proton-rich (Y e =t 0.52 - 0.54) for more than 10 seconds. 




(b, 



01 23456789 10 
Time After Bounce [s] 



X 






dyn 




X 




P 




X 




QW 











01 23456789 10 
Time After Bounce [s] 




(b) 



01 23456789 10 
Time After Bounce [s] 



T 






dyn 




X 




P 








T QW 









0123456789 
Time After Bounce [s] 



Figure 19: Comparing the mass accretion rate from our simulations and the approximation from Qian and Woosley 
( 1996) in the graphs (a) and different dynamic timescale in the graphs (b) with respect to time after bounce for the 10.8 
(left panel) and 18 M Q (right panel) progenitor models. 



8 Summary and Outlook 

For the first time, spherically symmetric core collapse supernova simulations based on general relativistic radiation 
hydrodynamics and three flavor Boltzmann neutrino transport are performed consistently for more than 20 seconds. We 
follow the dynamical evolution of low and intermediate mass progenitors through the collapse, bounce, post bounce, 
explosion and neutrino driven wind phases. The explosions of Fe-core progenitors of 10.8 and 18 M are modeled 
using artificially enhanced opacities, while the explosion of the 8.8 M O-Ne-Mg-core was obtained using the standard 
opacities. For all models under investigation, we confirm the formation and illustrate the conditions for the appearance 
of the neutrino driven wind during the dynamical evolution after the explosions have been launched. For the O-Ne- 
Mg-core and the 10.8 M Fe-core progenitor models, the supersonic neutrino driven wind collides with the slower 



expanding explosion ejecta where due to the deceleration a reverse shock forms, known as the neutrino driven wind 
termination shock. We discuss the impact of the reverse shock on several properties of the ejecta and find general 
agreement with Arcones et al. ( 2007| l. 

The comparison with approximate and static steady-state as well as parametrized dynamic wind models leads to a 
discrepancy in the obtained physical properties of the neutrino driven wind. Although the evolution of the hydrody- 
namic variables are in general agreement, we find smaller neutrino luminosities and a different behavior of the mean 
neutrino energies. In particular, the differences between the neutrino and antineutrino luminosities and mean neutrino 
energies are smaller. These differences reduce with time as the PNS contracts and result in a generally proton-rich 
neutrino driven wind over more than 10 seconds. Hence, the suggestion that the physical conditions in the neutrino 
driven wind could be favorable for the nucleosynthesis of heavy elements via the r-process could not be confirmed. For 
a more accurate determination of the yields of the neutrino driven wind, a detailed nucleosynthesis analysis taking the 
r-, p- and vp-processes into account is required. In order to further support the robustness, improvements of the input 
physics such as weak magnetism and nucleon-nucleon recoil (following e.g. Horowitz ( 2002| l), taking the presence of 
light and heavy clusters of nuclei into account as well as different EoSs with respect to the different PNS contraction 
behavior, should be considered. These may have a strong influence on the properties of the neutrino spectra at the 
neutrinospheres and may therefore modify some of the the results found in the present study of the neutrino driven 
wind. 

Our simulations are carried out until the neutrino driven wind settles down to a quasi-stationary state leading to 
the initial and neutrino dominated PNS cooling phase. There, the simulations have to be stopped because important 
neutrino cooling processes like the direct and modified URCA processes are not taken into account yet. However, a 
smooth connection to isolated neutron or protoneutron star cooling studies comes into reach for future work ( |Henderson| 
andPage|P007) ). 
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